Effect of in-plane magnetic field on magnetic phase transitions in u = 2 bilayer 

quantum Hall systems 
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By using the effective bosonic spin theory, which is recently proposed by Demler and Das Sarma 
[ Phys. Rev. Lett. 82, 3895 (1999) ], we analyze the effect of an external in-plane magnetic field on 
the magnetic phase transitions of the bilayer quantum Hall system at filling factor u — 2. It is found 
that the quantum phase diagram is modified by the in-plane magnetic field. Therefore, quantum 
phase transitions can be induced simply by tilting the magnetic field. The general behavior of 
the critical tilted angle for different layer separations and interlayer tunneling amplitudes is shown. 
We find that the critical tilted angles being calculated agree very well with the reported values. 
Moreover, a universal critical exponent for the transition from the canted antiferromagnetic phase 
to the ferromagnetic phase is found to be equal to 1/2 within the present effective theory. 

73.40. Hm, 73.20.Dx, 75.30.Kz 



Based on a microscopic Hartree-Fock theory and a long 
wavelength field theory,EJ recent theoretical works predict 
that, in the iy=2 bilayer quantum Hall (QH) system under 
quite general experimental conditions, there can be three 
qualitatively different quantum magnetic phases: the 
fully spin polarized ferromagnetic phase (F), the para- 
magnetic symmetric or spin singlet (S) phase, and the 
intermediate canted antiferromagnetic (C) phase. There 
have been some encouraging experimental evidence for 
the existence of the novel |C phase through inelastic list^ 
scattering spectroscopyj3u transport measurementsp'El 
and capacitance spectroscopy.lj However, it is difficult 
to make a precise experimental verification of the pro- 
posed quantum phase transitions, because a given sam- 
ple (with fixed values of the system parameters such as 
well widths, separations, etc.) is always at a fixed point 
in the quantum phase diagram calculated in Refs. [ 
Hence some exjtca, parameters (for exampJe, the in-plane 
magnetic fieldulj'tl or the bias voltagffltl'lD) are necessary 
in order to allow experimentally a continuous tuning of 
the ly — 2 QH system through the phase boundaries. 

In this paper, we present a theoretical investigation of 
the influence of an external in-plane magnetic field on 
the quantum phase diagram of the bilayer system. Al- 
though the topology of the phase diagram given by the 
Hartree-Fock theory should be correct, the Hartree-Epck 
theory overestimates the stability of the C phase. EJ In 
order to go beyond the Hartree-Fock theory, an effective 
hard-care boson theory was proposed by Demler and Das 
Sarma.B Later, it is pointed out that the hard-core bosons 
introduced in Ref. [ ^ form an exact Schwinger boson 
representation of an effective spin-1/2 system, and thus 
theoretical techniques developed for spin, models can be 
readily applied to the present problem Jj Hence we ex- 
tend this effective model to the case that a parallel mag- 
netic field is applied, and study its effect on the magnetic; 
phase transitions of the bilayer QH system a,t v — 2X3 
It is found that the system may make multiple transi- 



tions from the S phase to the C phase, and finally to 
the F phase as the parallel magnetic field increases. The 
general behavior of the critical tilted angle for different 
layer separations and interlayer tunneling amplitudes is 
shown, and we find that the critical tilted angles for the 
transition from the C phase to the F phase coincide quan- 
titatively with the experimental valuesJj'u Moreover, it is 
found that the critical exponent of this transition has a 
universal value, which is equal to 1/2 within the present 
model. 

In this work we assume that the electron-electron inter- 
action and Zeeman energies are much smaller than the 
Landau level splitting due to the strong perpendicular 
magnetic field B±z, and we therefore restrict the orbital 
Hilbert space to the lowest Landau level of either spin. 
When a parallel magnetic field along the y direction is 
added, with the choice of the gauge A = (0, B±x, — i3||a;), 
the Hamiltonian ofj-the double-layer system has the form 
H = Ho+Hi withtil 
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where c\ f. ^ creates an electron in the a-th layer 
(a = 1,2) with spatial wavefunction Uk{x,y) — 
^-iky^-{x^ki^f/2i' / ^^1/211 and spin a/2 (a = ±1) in 
the direction of the total magnetic field B = i?||y -f- B±z. 
[L is the length of system along the y direction and 
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I = y/hc/eB± is the niagne|tif; length.) Here the pa- 
rameters in Hq are defined ast£l 

AsAs = AsAsexp(-Q2;2/4)^ (3) 
A, = A,^l + (S||/Bi)2, (4) 

where Asas is the tunneling-induced symmetric- 
antisymmetric energy separation, and A^ is the Zeeman 
energy in the absence of -B|| . We notice that _B|| induces 
an Aharonov-Bohm phase factor exp{±iQkl'^) depending 
on the sense of interlayer tunneling with Q = B\^d/ B±J? . 

The matrix elements of the intralayer Coulomb inter- 
action are 

= ^«^(g)<5p,,,„e-^^'^/2e^^^^^'^ (5) 
q 

and the matrix elements of the interlayer Coulomb inter- 
action are 

Ve{puP2) = V'^ipi,P2) = V^\puP2) 

= ^^;s(g)<5p,,,„e-^^'^/2e'^^^^'^ (6) 
q 

Here VA{q) = {2t: / eq)FA{q,b) and VE{q) = 
VA{q)FE{q,b)e~'^'^ are the Fourier transforms of the in- 
tralayer and the interlayer Coulomb interaction poten- 
tials, respectively, e is the dielectric constant of the 
system, and d is the interlayer separation. In order 
to make a detailed comparison between the experimen- 
tal data and the theoretical predictions, we have also 
included finite-well-thickness correction by introducing 
the form factor FA{q,b) {FE{q,b)) in the intralayer (in- 
terlayer) Coulomb matrix elements, where FA(q,b) = 
2/bq - 2(1 - e-«'')/62g^ FE{q,b) = Asinh^ {qb/ 2) / b^ , 
and b is the width of a quantum well .113 

Due to the Aharonov-Bohm phase factor in the tun- 
neling process, if the electrons of the bilayer system are 
in the symmetric state for = 0, the inclusion of an 
in-plane field twists the original interlayer phase coher- 
ence aiid,|Pesults in an increase in the interlayer Coulomb 
energy.t^lEJ This effect can be made explicit, if one ab- 
sorbs the Aharonov-Bohm phase factor in Hq and makes 
the matrix elements of Hq real by redefining ci^k.a — 
exp{iQkl'^/2)ci^k,a and C2,fc,a- = exp(-z(5A;P/2)c2,fc,a-, 
which can be considered as a kind of pseudospin rotation. 
In terms of the new operators Ca,k,(n the matrix elements 
of the intralayer and the interlayer Coulomb interactions 
become 

t/"(pi,p2) = V^Hpi,P2) = Va{pi,P2), 

V'\p,,P2) = {V^\p^,P2)y = VE{pi,P2)e''^^''" ■ (7) 

Notice that the matrix elements of the interlayer 
Coulomb interaction become functions of Conse- 
quently, in term of the Ca^k,cT operators, the microscopic 



Hamiltonian at B\\ ^ can be considered as that at 
B\\ — with the modified matrix elements AsaSj A^, 
and V"'^"''^ . From these discussions, one can readily ex- 
tend the effective hard-core boson theory to the B\\ ^ 
case. We first j^ive a brief review of the model for the 
= caseJj'El and then provide the explicit relations 
between the matrix elements of the microscopic Hamil- 
tonian and the parameters in the effective hard-core bo- 
son theory, such that the effect of an external in-plane 
magnetic field can be easily incorporated. 

The Hamiltonian of a simple bilayer lattice model to 
describe the physics of the bilayer v = 2 QH system at 
Bjl ~ may be written ascl 

T, AsAs t , t ^ 

i 

i 

+ ^Y.i^nu-lf + {n2.^lf) 

i 

-J^(Si,;Si, +82,82,) (8) 

where i is the in-plane site (intra-Landau-level) index, 
and a is the spin index. S^^ = Y^^p Aia 

and nii — c^ia^iia are spin and charge operators for 
layer 1, with analogous definitions for layer 2. The effec- 
tive Heisenberg coupling J and the local charging energy 
tc of this model can be estimated as follows. In order to 
have the sarrua-magnon spectrum in the F phase, one must 
impose thatEj Jo? /2 = Al^ with the lattice constant 
a = V2^l and A = (1/4) WA(fc)|fcpZ2 exp(-fc2p/2). 
Thus 

J=^Y.^A{k)m'e->^'^'/'. (9) 

k 

On the other hand, Cc can be estimated from Hi under 
the Hartree-Fock approximation for the system in the F 
phase, which is given by 

ec = J2VA{k,0)-Y,VE{k,0). (10) 

k k 

Under the simplifications in which the total charge 
fluctuations are left out and only the lowest two energy 
states for a given Landau orbital are kept, the effective bi- 
layer lattice ropdel can be further reduced to a hard-core 
boson theory,El and it leads to a phase diagram which is 
more precise than that given by the Hartree-Fock theorw 
and is actually exact within the reduced Hamiltonian.EI 
The phase boundary separating the F and the C phases 
and that separating the C and the S phases can be writ- 
ten as 

A;, = - J(l - sin26'), (11) 
= - J(H-sin26l), (12) 
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where the parameters Ey and 9 are 



Ey — 



e = tan^^ 



^Ias + (^c/2)2, 
ec/2 



iSAS 



2 

SAS 



(13) 
(14) 



Now we turn to the case of the presence of an in- 
plane magnetic field. As mentioned before, when ^ 
0, one needs to replace Asas and by Asas and 
A2, respectively. Moreover, the matrix elements of 
the interlayer Coulomb interaction become Ve{pi,P2) — 
Ve{pi,P2) Gxp{±ipiQP). Therefore, the local charging 
energy becomes 



2 roc 

ec = ^ dxe-^"/^FA{x/l,b) 
Jo 



2 roo 



^ dxe-"" ^^e-'"^^^FE{x/l,b)Jo{Qxl) (15) 



el 



{ Jo{x) is the-|Bessel function), which is an increasing func- 
tion of 

Substituting these modified parameters into Eqs. (pj])- 
(p^), the phase boundaries at nonzero are obtained. 
The quantum phase diagram for several in-plane mag- 
netic fields is shown in Fig. |^. (In the foUowings, the 
length and the energy units are chosen to be the mag- 
netic length / and the intralayer Coulomb energy /el.) 
It is clear that the region of the F (S) phase is expanded 
(shrunk) as B\\ increases. The reason that the F phase 
is enhanced and the S phase is suppressed comes from 
the fact that the system at B\\ 7^ has a weaker Asas 
and a stronger A^ compared to the system at B\\ — 0. 
From this quantum phase diagram, one finds that a par- 
allel magnetic field can be a useful control parameter to 
tune the system through the phase boundaries. For ex- 
ample, by applying _B|| , the system can undergo a phase 
transition from the C phase to the F phase (or from the 
S phase to the C phase). Indeed, such transitions are 
observed in recent tilted-field experimcnts.H'E[Q The sam- 
ples in Refs. [ |^ and [ ^ are initially in the C phase near 
the F-C phase boundary in the absence of B\\. The bi- 
layer system transits to the F phase when the tilted angle 
Q = tan~^(_B||/_Bj_) reaches a certain critical value 9c. 
Thus the present theory can provide a qualitative expla- 
nation of the observed transitions in these experiments. 

Moreover, we are able to compare the calculated &c 
with these experimental results quantitatively. Since it 
is relatively easy to vary AgAs and d in fabrication, we 
focus our attention on the dependence of Qc on them. 
Our theoretical prediction of the dependence of Qc on 
Asas (while the other parameters are fixed) is shown in 
Fig. ||. We find that Qc is a monotonic increasing func- 
tion of Asas- It is expected because the system with a 
larger Asas locates at the phase diagram with a farther 
distance to the F-C phase boundary, and a larger bend- 
ing of the F-C phase boundary caused by a bigger is 



necessary to make a transition to the F phase. Adopt- 
ing the sample parameters reported in Refs.[ ^ and [ |), 
we obtain the theoretical values of Qc — 33.5° and 48°, 
respectively (see Fig. |^), which agree very well with the 
measured values: Qc — 37° in Ref.[ ||] and Qc — 50° 
in Ref. [ ^. The close agreement gives us confidence on 
the accuracy of the phase diagram calculated by using 
the effective bosonic spin theory. We also find that, near 
the transition point Asas = ^SAS' both curves in Fig. 
(and other curves for several sets of system parameters^ 
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fit very well with the function Oc — (Asas ~ A. 
where a = 1/2. Indeed, by using series expansion for 

small Qc and Asas ~ ^sas' this functional form can be 
derived from the formula of the F-C phase boundary. 
This result means that a is a universal critical exponent 
at the C-F phase transition, and its value (a = 1/2) may 
reflect the mean-field character of the present approach. 

The dependence of the critical tilted angle Qc on layer 
separation d is shown in Fig. ^ for several different sets of 
sample parameters. It can be seen that the critical tilted 
angle drops slowly at large d. In fact, for some cases it 
seems that Qc may not drop to zero even when the layer 
separation approaches infinity, which means that the C 
phase is stabilized in the single-layer limit. This is an 
artifact due to fixing the interlayer tunneling amplitude 
Asas for a given curve. If the dependence of Asas(c?) 
can be known and incorporated in the calculation, the 
critical tilted angles should eventually drop to zero since 
the C phase is expected to disappear when the two layers 
are completely decoupled. Therefore, the increase of d in 
Fig. ^ should not be interpreted simply as a process of 
moving two layers apart. Rigorously speaking, a single 
curve in Fig. ^ is merely meant for samples with different 
values of d but the same value of Asas- 

In conclusion, we have shown that the effective bosonic 
spin theory with no adjustable parameters can describe 
the experimental observations for the quantum phase 
transitions with quantitative accuracy. Moreover, we find 
that the critical exponent of the C-F phase transition has 
a universal value a = 1/2 within the present effective 
model. It would be quite interesting to see whether fu- 
ture experiments support our prediction or not. Finally, 
the results presented here can be useful guidelines for ex- 
perimentalists to design their samples for observing the 
magnetic phase transitions by tilting magnetic fields. 
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FIG. 1. 1^=2 bilayer phase diagrams in the bosonic 
spin theory for different in-plane magnetic fields . Here 
Az means the Zeeman energy caused by the total magnetic 
field. Continuous, dotted, and dashed lines correspond to 
B^\/B± = 0, l/\/3, and 1, respectively. The width of the 
electron layer is 1.0 and the interlayer separation is 1.45. 
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FIG. 2. Critical tilted angles for the C-F phase tran- 
sition are plotted as a function of Asas. Circles are ob- 
tained using the parameters for the sample in Ref.[ ^, where 
{Az,d,b) = (0.008,1.45,1.0); squares are for the sample in 
Ref.[ |, where (Az,d,6) = (0.00687,1.08,0.94). To com- 
pare the critical tilted angles with experimental values, choose 
Asas = 0.10 for the first sample and Asas = 0.117 for the 
second sample. 
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FIG. 3. Critical tilted angles for the C-F phase transi- 
tion are plotted as a function of layer separation d. All three 
curves have the same A2(= 0.008) and b (=1.0) The values 
of Asas for circles, squares, and diamonds are 0.1, 0.09, and 
0.08 respectively. 
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